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The complex arrangements of atoms near grain boundaries are difficult to understand 
theoretically. We propose a phenomenological (Ginzburg-Landau-like) description of 
crystalline phases based on symmetries and fairly general stability arguments. This 
method allows a very detailed description of defects at the lattice scale with virtually no 
tunning parameters, unlike usual phase-field methods. The model equations are directly 
inspired from those used in a very different physical context, namely, the formation of 
periodic patterns in systems out-of-equilibrium (e.g. Rayleigh-Bcnard convection, Turing 
patterns). We apply the formalism to the study of symmetric tilt boundaries. Our results 
are in quantitative agreement with those predicted by a recent crystallographic theory 
of grain boundaries based on a geometrical quasicrystal-like construction. These results 
suggest that frustration and competition effects near defects in crystalline arrangements 
have some universal features, of interest in solids or other periodic phases. 

Keywords: Grain boundaries; Pattern formation; Geometrical frustration; Quasicrystals; 
Faceting. 



1. Introduction 

The structure and dynamics of grain boundaries in materials are still little 
understood. 1,2 The possible relationships between the mesoscopic features of grain 
boundaries and their inner microscopic structure at the atomic scale remain un- 
clear. The energies of some interfaces can be accurately described with continuous 
theories and topological arguments. But the study of other quantities (like mo- 
bilities), as well as short scale phenomena such as faceting, may require a deeper 
understanding of the microscopic structure, which directly reflects the discrete and 
anisotropic character of the lattice. Near a boundary, the crystal structure usually 
has no simple periodic behavior and is therefore difficult to characterize. Any quan- 
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titative and sufficiently simple theory allowing a systematic characterization of the 
fine structure of grain boundaries at zero temperature could open the way toward 
a general description of these objects. Following such an approach requires to go 
beyond usual continuous media approximations. 

In this paper, we discuss a way of elucidating grain boundary microscopic struc- 
tures within a theoretical framework that explicitly takes into account spatial vari- 
ations at the atomic length scale. The models considered here have a very small 
number of parameters. They derive from very general energy and symmetry argu- 
ments. Although the nonlinear approach proposed below is fairly uncommon for 
crystals, it is largely used in a (apparently) completely different field of Physics, 
namely, in pattern formation problems. Many physical systems driven out of equi- 
librium spontaneously form periodic structures, some of which are characterized by 
a crystal-like order. 3 Some well-known examples are the cellular structures formed 
in Rayleigh-Benard convection of heated fluid layers, or Turing patterns in reaction- 
diffusion processes. Pattern formation is ubiquitous in nature and has universal fea- 
tures. Some of these features may be of relevance for solid crystals as well. Striking 
analogies between pattern forming systems and dissipative processes taking place 
in solid crystals have been recently pointed out. 4 ' 5 In Section 2, we briefly justify 
why some simple nonlinear models of pattern formation may be useful to study 
defects in crystals. We show in Sec. 3 that these models are very efficient to cap- 
ture elementary competition and geometrical exclusion effects that probably take 
place in genuine grain boundaries. In particular, they are able to predict interface 
faceting, a feature still little understood from continuous theories. 6 The results ob- 
tained arc in very good agreement with those given independently by a theory of 
grain boundaries developed recently, which is based on geometrical arguments and 
inspired from quasicrystals theory. 7 

2. Basics of pattern formation 

We briefly recall some concepts of pattern formation, that can be found in well doc- 
umented reviews or textbooks. 3 ' 8 ' 9 Consider a spatially extended system described 
by a phenomcnological, dimcnsionlcss local order parameter xfj(x,t) depending a 
priori on space and time variables. In convection problems, for instance, tp(x, t) is 
related to the vertical velocity at the mid-plane of the convective cell. A uniform 
phase is characterized by the trivial solution ip =cst(= 0). When a non-trivial sta- 
tionary periodic pattern sets up in the system, e.g. above the onset of convection, 
ip can be written in a first approximation as a sum of plane waves: 

N 
n=l 

where N is the number of wavevectors (of same modulus \k n \ = k = 2-k/Xq) 
characterizing the base pattern (in two spatial dimensions, N = 3 for hexagonal 
symmetry and TV = 2 for square symmetry). A perfect crystal is observed when 
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the amplitudes A n are all equal to a constant A . The Figures 2 and 3 further 
show examples of the field ip in gray scale in 2D (dark regions have tp(x) > 0, 
bright ones ip < 0). We now look for a simple phenomenological isotropic partial 
differential equation satisfied by ip. Keeping in mind the analogy with crystals, in 
addition to the periodicity Ao, one wishes an equation that imposes a symmetry on 
the stationary solution (for instance hexagonal, TV = 3). A linear wave equation 
like (A + kfyip = is not acceptable, as a superposition of any N > 1 waves 
would be solution. Instead, we need to resort to a nonlinear equation, introducing a 
coupling between the different amplitudes A n in Eq.(l). We consider here the well- 
known Swift- Hohcnberg model (first introduced in the context of Rayleigh-Benard 
convection), which also involve the time variable t: 3 



where e <C 1 is a dimcnsionlcss control parameter, the reduced Rayleigh number 
in convection problems. If e < in (2), ip = is the only stationary solution 
(dip/dt = 0) that is stable. The uniform solution ip = becomes unstable for 
e > 0. In that case, new solutions with broken symmetries appear, in the form of a 
periodic pattern characterized by a layer spacing Ao = 2-7r/fco. The order parameter 
ip can be approximated by a sum of the form (1) provided that e 1. Nonlinear 
terms are responsible for mode combinations: given an arbitrary set of wavevectors 
{fe„}, one can derive equations for their amplitudes A n , these amplitudes being now 
coupled to each other. The A„'s happen to have stable non-zero solutions only in 
a very small number of cases, that generally correspond to symmetric, crystal-like 
arrangements of the {k n }. For instance, equation (2) with g 2 = leads to stripes 
formation (N = 1). On the other hand, if g 2 ^ and < e < 4g|/3 only hexagonal 
solutions are stable. 9 

Let us now point out one of the crucial properties of this type of model: the 
evolution equation (2) has a potential structure. One can associate an "energy" , or 
Liapunov functional F, to the system, such that Eq.(2) reads dip/dt = —SF/Sip 
(the right-hand-side representing a functional derivative) . F is given in that case by 



It is easy to show from Eq.(2) that dF/dt < 0: the quantity F always decreases with 
time. The dynamics is dissipative and converges toward stationary stable solutions 
that correspond to local minima of F in the configuration space. The absolute mini- 
mum of F corresponds a perfectly regular pattern of the form (1), with wavenumber 
k{) (a "perfect crystal"). Consider now an immobile grain boundary separating two 
differently oriented domains that arc otherwise in the ground state: near the bound- 
ary, it is impossible to preserve the regular crystal structure imposed in the bulk. 
As a consequence, such a configuration must have a higher F. 

One can now see the pattern formation process given by (2) from an alternate 
view-point, very useful for applications to materials: namely, as an energy minimiza- 



# = ^-i(fc 2 + A) 2 V + ^ 2 -^ 3 , 



(2) 



(3) 
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tion problem within some constraints. Stationary grain boundaries can be regarded 
as optimal solutions of a frustration problem of two competing lattices of different 
orientations. The free-energy functional imposes a fixed periodicity on the field tp in 
the bulk, minimizing the term ?/;(£; 2 + A) 2 -0 in (3), as well as a given bulk symmetry, 
which is fixed by mode combination through the higher order terms. The equations 
for the amplitudes A n show that the grain boundary is actually an interference 
pattern between the two crystals, with amplitudes not constant in space and not 
equal to each other. 

Different choices of nonlinear terms can lead to different crystal symmetries in 
the bulk. Adding other cubic or higher order terms in the functional (3) generates 
new nonlinear terms in Eq.(2), modifying the relative coupling constants between 
the different amplitudes A n in the solution (1). As a consequence, a crystal formerly 
unstable may become stable. For instance, a square crystal lattice (N — 2) can be 
obtained by choosing: 10 F = J dx[-{e/2)i}j 2 + l/(2fc^)-0(fcg + A) 2 ^ - ff 2 ^ 3 / 3 + 
/?V 4 /4+ (7/4)^ 2 A 2 (?A 2 )], with 13 < and e > 16^(3/3+ 13 7 )/(3/3 + 4 7 ) 2 . This 
free energy gives coefficients that depend on the angle between the wavevectors in 
the amplitude equations, such that they stabilize solutions with perpendicular wave 
numbers. 

The energy (3) involves a Taylor series expansion of an order parameter, as usual 
in Ginzburg-Landau descriptions of phase transitions. This simple form suggests 
that geometrical frustration in crystalline phases with isotropic interactions is likely 
to have some universal features. Actually, numerical solutions of equation (2), or 
slight variants of it, are able to reproduce known features of crystals. The behavior 
of the energy F of grain boundaries as a function of misorientation angle follows 
a Read-Shockley law. 4 In addition, short-range pinning forces act on defects, a 
phenomenon analogous to the Peierls stress in crystals. 5 

In the following Section, we emphasize the close relationships that exist between 
this pattern formation approach of grain boundaries and some space filling prob- 
lems under constraints, that have been solved in Ref. 7 by the use of quasi-crystal 
methods. We briefly recall the main aspects of this second theory below. 

3. Structures of grain boundaries 
3.1. A geometrical construction 

The main distinctive feature of this theory is that it considers interfaces and qua- 
sicrystals as a region in space where two or more interpenetrating crystal lattices 
compete for space. The final atomic positions are then decided by a modified version 
of the strip-projection proposed by Katz and Duneau 11,12 to study quasicrystals. 
Perhaps the most important attribute of this approach is that it allows interfaces 
and quasicrystals to be described by the same set of equations, thus rendering them 
formally equivalent. 

Given a completely arbitrary set of lattices and relative orientations the formal- 
ism produces ideal (minimum local strain) structures which are expected to play for 
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Fig. 1. Schematic representation of 6D space showing the orthogonal 3d subspaces Vi, Vs, i?" , E 1 - 
as ID lines. The crystal lattices Li and L2 are respectively contained in V\,Va. The interface is 
defined as the set of points cc" that are projections of the hyperlattice points located inside a small 
region around fill. This region is called the strip and is delimited by the two thin lines (see text). 

grain boundaries the same role than the perfect crystal concept plays for crystals, 
i.e., they define defect free, reference structures analogous to Bravais lattices. Since 
strain minimization is a physical consideration, the formalism is endowed with a 
physical basis in spite of its geometrical formulation. 

Full details of this method can be found in Ref. 7. Briefly, the crystal lattices are 
embedded in a higher dimensional space where they define an hyper lattice whose 
symmetry depends on that of the lattices and their relative orientation. This hyper 
lattice is then projected back to physical space to produce the actual interfacial 
structure or quasicrystal. 

The higher dimensional 61? space is schematically represented in Figure 1. It 
is divided in two sets of orthogonal 3D subspaces: (Vi,V2) and (E\E ± ). Vi and 
V2 contain the lattices L\ and L 2 of the two crystals, and E», E 1 - contain respec- 
tively, the actual interfacial structure and the strain field across the interface. E L 
is actually the Bollmann's displacement space or b-space. 1 

Only those hyperpoints contained within a small region around E" (see Figure 
1) called the strip 11,12 are projected. The size of the strip is chosen to include pairs 
of atoms that occupy incompatible positions, these atoms are then replaced in E^ 
by a single atom at their average position. Given two interpenetrating lattices in 
physical space, the method constructs an ideal, best fit, minimum strain lattice as 
the set of points of E" . By adopting the average position, atoms at the interface act 
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as a strain buffer between the crystals on each side of the interface. The component 
in E 1 - is a measure of the local strain or frustration between two nearly coincident 
(overlapping) atoms from each lattice. 

This method allows the construction of a complete bicrystal with the parent 
crystals separated by a unique, fully 3D structure called the interfacial layer (IL), 
arising from the projection of hyperpoints within the strip. The role of the IL, whose 
structure is different from that of either crystal is to resolve the positional incompat- 
ibilities that occur at the interface. Since grain boundaries are usually 2D, a suitable 
IL plane (or surface) must be selected in order to model any particular boundary 
plane, of given relative orientation to the parent crystals. A bicrystal can then be 
built by filling the space on cither side of the boundary with crystal sites projected 
from (Vi) and (V2). Note this approach is unique in that the grain boundary always 
contains atomic (IL) sites, which makes it susceptible of crystallographic analysis. 

We recall now the definition of the O-lattice, a concept that will be useful in 
the next Section. Given two crystal lattices, L\, L 2 and a transformation A, such 
that L 2 = Aii, in the absence of translations, the O-lattice is defined as the set of 
points that satisfy 

(I-A- 1 )- 1 ^ 1 ) =z(°) (4) 

where I is the identity matrix, and G L\. Geometrically, the O-lattice corre- 
sponds to points in space (not necessarily lattice points) of minimum strain i.e., 
minimum misfit between L\ and L 2 - Another interesting property of the O-lattice 
is that if ZW* and (A/W) with i = 1,2,3 are the base vectors of the reciprocal 
lattices of L\ and L 2 respectively, then 

(/W*-(AZW)*).x(°)=n (5) 

with integer n. This equation establishes a relationship between the O-lattice and 
the Moire pattern formed by the superposition of L\ and L 2 . 

3.2. Results 

In this section we study symmetric tilt grain boundaries in two-dimensional (2D) 
crystals. The initial condition ip(x,t = 0) in Eq.(2), or its analogue for square 
patterns, is made of two crystals rotated by 9 m i S , the interface being directed along 
the vertical axis (see Ref. 5 for more details). We introduce an angle 9 that either 
corresponds to the misoricntation angle m i S or to 90° — 9 m is (for square patterns) , 
depending on the orientation of the symmetric boundary plane. Note that < 
9 mis < 45° and < 9 < 90°. For each rotation 9 m i S , there are two possible boundary 
orientations (with different structure) that lead to a symmetric interface. The angle 
9 represents the angle formed between two atomic planes that are symmetric with 
respect to the boundary. 

We analyze the stationary (minimum energy) fields obtained at large times by 
solving Eq.(2) numerically, and compare their structures with those given by the 
method described in Sec. 3.1. 
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Fig. 2. Singular tilt boundary with 9 = 28.072° (S17). Right: if) in gray scale given by the 
nonlinear model. Left: from the construction of Sec. 3.1. See text for legend. 

Figure 2 displays (on the right) a boundary obtained by solving the Swift- 
Hohenberg equation for square patterns, with an angle 9 — 2 tan _1 (l/n), where n 
is an integer (4 in this case). Such interfaces are singular 7 and characterized by a 
simple repeating pattern, as emphasized by the bright circles. 

In the geometrical method, any grain boundary structure is a 2D manifold (a 
plane, usually) included in the 3D IL. The most likely manifolds are a priori the ones 
that run along directions that maximize the density of IL sites. Such dense planes 
must clearly intersect perpendicular IL planes along directions that have also a high 
density of sites. In the particular case of a tilt boundary configuration, the planes of 
the IL normal to the tilt boundary constitute the model's prediction for the so called 
twist boundary. Hence, in the geometrical approach, twist boundaries do provide 
information about the possible directions that a perpendicular (tilt) boundary may 
run along. We will hence focus on the structure of twist planes below. 

The left part of Figure 2 shows the twist interface corresponding to the example 
above, obtained using the approach of Sec. 3.1 for a cubic lattice. It is a plane normal 
to the (001) direction in the cubic (100)/ (001) twist interface. As indicated above, 
this plane is a map of all the possible transverse conformations (or trajectories) of a 
tilt grain boundary between two grains of given misorientation angle. A particular 
trajectory is indicated in gray. Large disks represent atomic positions, or interfacial 
sites. The short and long arrows indicate the O-lattice (OL) and coincidence sites 
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lattice (CSL) base vectors, respectively. Note that the interface is composed of 
identical domains centered at O-points (small dots), that are points of zero strain. 
Atomic domains are separated by a network of primary screw dislocations. 1 The 
intersection of these dislocations with the boundary plane is highlighted by a square 
grid in the figure. 

The predicted atomic positions (grayed dots) match nearly exactly the (black) 
regions where ip is maximum in the nonlinear model. The short horizontal lines 
on the figure mark the positions of dislocations. These are edge dislocations, with 
Burger's vectors normal to the boundary line. A similar remarkable agreement was 
obtained in every case for a variety of angles 9 chosen arbitrary, in square and 
hexagonal symmetries. Non-singular interfaces are characterized by more compli- 
cated interfacial patterns (not shown). 

Since elastic strain is minimized near O-points, it is expected that tilt interfaces 
of low angles 9 tend to run along directions that follow the O-lattice, preserving this 
way the structure of the domains around the O-points. The O-lattice has the same 
symmetry as the crystal (see Figure 2, left), but a distinct period Ao, given by Ao = 
A/ [2 sin(0/2)] independently of the symmetry 7 (A is the period of the crystal; Ao > A 
at low angles). Therefore, the grain boundaries separating two symmetric grains 
along the vertical axis should be located at some discrete x-positions, separated from 
each other by Ao (Figure 2, left). On the other hand, a weakly nonlinear analysis 
of Eq.(2) allows the derivation of the law of motion of straight grain boundaries. 5 
It is found analytically that interfaces feel a pinning potential, periodic along the 
x-direction. Strikingly, the nonlinear model agrees perfectly with the geometrical 
construction: The periodicity of the pinning potential is A/[2sin(0/2)] and coincides 
with Ao- In the nonlinear model, the derivation of this formula follows from a 
solubility condition that imposes particular combinations of the modes k n . The 
result is formally equivalent to the the formula (5) above. The lines going through 
points of the O-lattice thus correspond to minima of a pinning potential in the 
nonlinear model. 

Figure 3 (on the right) represents a configuration obtained at time t = 2000 
from the same nonlinear model, with a high angle 9 = 73.74°. The initially straight 
interface quickly destabilizes into a broken-symmetry phase, characterized by facets 
of short period. This faceting instability is not observed for low 9 boundaries, for 
instance that of 9 = 90 — 73.74, the other symmetric boundary orientation in this 
bicrystal. 

The faceted structure is also visible from the geometrical construction of 
(left). In this configuration, a straight tilt grain boundary would run vertically, i.e. 
along the diagonal of the domains of best fit of size 3 or 4, where the distance 
between O-points is larger. In order to keep the structure of these domains any- 
way, the interface becomes faceted instead. It runs along the core lines of the screw 
dislocation network. In this case, the interface is locally aligned with the Burgers 
vectors of the dislocations that compose it. Along this line, a repeating structure of 
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Fig. 3. Facet formation. Right: Nonlinear model with at 9 = 73.74° (S25). Left: Geometrical 
construction. 

4 atoms facing 3 atoms of the other grain appears, as correctly reproduced by the 
nonlinear method. Note that the geometrical map of Figure 3 (left) also contains 
the structure of symmetric tilt boundaries with the angle 9 = 90 — 73.74 = 16.26°. 
The direction of these interfaces is now inclined of 45° with respect to the vertical 
direction in Figure 3 (left): the most likely trajectories are straight and run along a 
line of maximum O-point density (or pinning potential minimum) , as confirmed by 
nonlinear calculations (not shown) . If an interface with such an orientation destabi- 
lized, it would find it difficult to develop the same "Chinese hat" facets as observed 
for the other orientation 9 = 73.74°. Hence, low 9 grain boundaries are probably 
less likely to facet. In conclusion, when the average orientation of a boundary is not 
parallel to the O-lattice, the boundary may become faceted. 

Preliminary numerical results on the nonlinear model for square patterns show 
that the faceting instability does not occur for misorientation angles lower than a 
characteristic value 6 C w 60°. This value can be interpreted qualitatively as follows: 
For 9 > 6 C , the periodicity of the pinning potential A/[2sin(0/2)] becomes similar 
to that of the crystal making it easier for the interface to deform under stress. 

The configuration of Figure 3 (right) has a relatively high elastic stress. After 
sufficiently long times, it starts to slowly evolve into more complicated patterns, 
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with longer facet lengths in order to relieve the initial stress. 
4. Conclusions 

We have presented a theory of crystals based on a local order parameter description 
and general stability and symmetry principles. The theory is well suited for the 
description of short scale phenomena, in particular the intricate spatial distribu- 
tions of atoms in grain boundaries. The results compare very well with those given 
independently by a generic geometrical theory of grain boundaries. The nonlinear 
nature of the model naturally incorporate some underlying topological structures, 
such as the O-lattice, which in this theory is interpreted as the set of boundary 
pinning sites. 
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